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Abstract 

Nested Archimedean copulas recently gained interest since they generalize the well- 
known class of Archimedean copulas to allow for partial asymmetry. Sampling 
algorithms and strategies have been well investigated for nested Archimedean copulas. 
However, for likelihood based inference it is important to have the density. The 
present work fills this gap. A general formula for the derivatives of the nodes and 
inner generators appearing in nested Archimedean copulas is developed. This leads 
to a tractable formula for the density of nested Archimedean copulas in arbitrary 
dimensions if the number of nesting levels is not too large. Various examples including 
famous Archimedean families and transformations of such are given. Furthermore, a 
numerically efficient way to evaluate the log-density is presented. 
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1 Introduction 

£N) There has recently been interest in multivariate hierarchical models, that is, models that 

are able to capture different dependencies between and within different groups of random 
variables. One such class of models is based on nested Archimedean copulas. A partially 
nested Archimedean copula C with two nesting levels and do child copulas (or sectors or 
groups), is given by 

O C(u) = Co(CiOi), . . . ,C do (u do )), u = («!, . . . ,u do ) T , (1) 

where do denotes the dimension of Co and each copula C s , s 6 {0, . . . , do}, is Archimedean 
with a completely monotone generator tb s , that is, 
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where 



t s (u s ) = ^1Ps 1 ( u sj) 



and ijj s : [0,oo] — > [0,1] is continuous, ip s (0) — 1) ips(oo) = Hindoo ip s it) = 0, and 
(— l) k t/ji k \t) > for all 6 No, t 6 (0,oo). The set of all completely monotone 
Archimedean generators is denoted by in what follows. The copula Co is referred to 
as root copula. Model ([I]) provides an intuitive hierarchical structure, since, for example, 
if U ~ C the pair (U s j, U s k) T (j ^ k) has joint copula C s whereas the pair (U r j, U s k) T 
(r ^ s) follows the root copula Co- One can therefore directly say what the bivariate 
margins are and theoretical results about measures of association, for example, directly 
apply. Furthermore, such a construction provides an explicit form for the copula itself, 
which is important, for example, in likelihood-based inference for censored data. More 
complicated nesting structures can be constructed; see Section [6j In general, a nested 
Archimedean copula is an Archimedean copula with arguments possibly replaced by other 
nested Archimedean copulas. 

For statistical applications it is desirable to be able to evaluate the density of a 
multivariate model, for example, for parameter estimation or when conditional copulas 
are required (for example, for goodness-of-fit testing via Rosenblatt's transform; see 



Genest et al. (20091). For Archimedean copulas, the density (if it exists) is theoretically 



trivial to write down; for one obtains 

C S (U S ) = 4 d) (ts(Us)) riCC'/K-), U ° G (°> 1 ) dS ' 

However, the appearing generator derivatives tpi are non-trivial to access theoretically 
and, even more, computationally. This issue has recently been solved for several well- 
known Archimedean copulas and transformations of such; see Hofert et al. ( 2012a[ ) and 



Hofert et al. (2012b). Our goal is to extend these results to the corresponding nested 
Archimedean copulas. Note that this is more challenging because differentiating is 
more complicated due to the inner derivatives that appear when applying the Chain 
Rule; in contrast to Archimedean copulas, these inner derivatives depend on variables 
with respect to which one has to differentiate again. Already in low dimensions the 
corresponding formulas for the density c become challenging to write down and, even 
more, to evaluate in a numerically stable way. 

After a brief introduction and overview about nested Archimedean copulas in Section [2j 
we tackle the problem of computing the density of ([T]) by first deriving a convenient form 
for the copula. This will allow us to compute the density; see Section |3j All necessary 
details for several well-known Archimedean families are provided in Section [4j Section 
[5] addresses numerical evaluation of the log-density. Section [6] presents the density for 
three-level nested Archimedean copulas and extensions to higher nesting levels are briefly 
addressed. For the reader's convenience, proofs are deferred to the appendix. 
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2 Nested Archimedean copulas 



The class of nested Archimedean copulas was first considered in Joe (1997, p. 87) in the 



three- and four-dimensional case and later by McNeil ( 2008 ) in the general d-dimensional 
case. McNeil (2008) and Hofert (2011a) derive an explicit stochastic representation 



for nested Archimedean copulas which allows for a fast sampling algorithm of nested 
Archimedean copulas similar to the Marshall-Olkin algorithm for Archimedean copulas; 



see Marshall and Olkin (1988) for the latter. Hofert (2011b) provides efficient sampling 



strategies for the most important ingredients to this algorithm, the random variables 
responsible for introducing hierarchical dependencies. An implementation for several 
well-known Archimedean families (and transformations of such) is provided by the R 



package copula; see Hofert and Machler (2011). 



Although nesting is possible in more complicated ways (see Section [6]) , in the following 
we focus on nested Archimedean copulas of Type ([TJ (with some child copulas possibly 
shrunk to single arguments of Co). By Bernstein's Theorem, each tp S is the Laplace— 
Stieltjes transform of a distribution function F on [0, oo) with -F(O) = 0. A sufficient 
condition under which ([I]) is indeed a proper copula is then that the nodes 



Qs 



ipQ 1 o ip s , s G {1, . . .,d }, 



have completely monotone first order derivatives; see McNeil (2008). Note that this 



sufficient nesting condition is indeed only sufficient but not necessary. For example, if 
ipo{t) = — log(l — (1 — e~ 6 °) exp(— t))/9o denotes the generator of a Frank copula and 
ipl(t) = (1 + i) -1 / 6 * 1 the generator of a Clayton copula, then C(u) = Cq(ui, C\(u2, U3)) 



-9 



)-l < 



is a valid (nested Archimedean) copula for all 9o,9\ such that 9q/(1 ■ 
although ipQi is not completely monotone for any parameters 9$,9\. 

Among the most widely used parametric Archimedean families are those of Ali-Mikhail- 



Haq, Clayton, Frank, Gumbel, and Joe; see Hofert et al. (2012b) for the corresponding 



generators, their derivatives, Laplace-Stieltjes inverses, and properties of the copula fami- 
lies. These one-parameter families can easily be extended to allow for more parameters, 
for example, via outer power transformations. For more details on this and other aspects 



of nested Archimedean copulas we refer to Hofert (2010) and the references therein. 



3 Inner generator derivatives and densities for two-level nested 
Archimedean copulas 

3.1 The basic idea 

Let C be a (f-dimensional nested Archimedean copula of Type ([I]) (with some child copulas 
possibly shrunk to single arguments of Co) and assume the sufficient nesting condition to 
hold; for the Ali-Mikhail-Haq, Clayton, Frank, Gumbel, and Joe families, this is fulfilled 
as long as all generators belong to the same family and #0 < 9 S , s G {1, . . . , do}. This 
condition implies that each copula C s , s G {1, . .. ,do}, is more concordant than Cq. 
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One of the main ingredients we need in the following is the function 

ipOs(t; v) = exp(-vipo s (t)) (3) 

which we refer to as inner generator. It is a proper generator in t for each v > 
as a composition of the completely monotone function exp(— v •) with ipQ S which has 
completely monotone derivative. With Fq = CS~ 1 [i(jq\, we obtain 

C(u) =C Q (Ci(u 1 ),...,C do (u rfo )) = exp^-voY^^Os(ts{us))jdFo(v ) 

poo d a 

= / Y[ijJo s (ts{u s );vo)dFo{v ) (4) 

J „ 1 



s=l 

By our assumption of having completely monotone generators, the density c of C exists 
and is given by 

Qd 

c(u) = o— C « • 

du doddo ■ ■ ■ du n 

Instead of differentiating ([!]) directly, the idea is now to use Representation Q. By 
differentiating under the integral sign, the density c allows for the representation 



roo d o d o d s 

c(u) = / Y[4 d ;\t s (u s );vo)dF (v ) ■ [] n^/K-) 

J ° s=l s=l j=l 



E 



do 



Y[4 d ;\t s (u s );V ) •niK^ 1 )'^)- (5) 

= 1 J 8=1 j=l 



do d s 



For the cost of one integral (which will be computed explicitly below!), one can therefore 
easily compute the density c (theoretically) as a Fo-mixture. This is especially advanta- 
geous in large dimensions as the complexity of the problem does not (again, theoretically) 
depend on the sizes of the child copulas too much, rather on the number of children. 
From Equation we identify the following key challenges: 

Challenge 1 Find the derivatives of the inner generators ipo s (t; vo); 
Challenge 2 Compute their product; 

Challenge 3 Integrate it with respect to the mixture distribution function Fq = [i/jq] . 



All three challenges will be solved in Section 3.3 with the help of the tools presented in 
the following section. 

3.2 The tools needed: Faa di Bruno's formula and Bell polynomials 

One formula which proves to be useful here, is the expression of the nth derivative of 



a composition of functions; see Craik (2005). Although this formula dates back to the 
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work of Arbogast (1800), it is named after the mathematician Faa di Bruno. For suitable 
functions / and g, Faa di Bruno's formula states that 



n 



k=i jev n!k \ J1, --- ,Jn - k+1 / i=i v L - 1 

where ( . n . ) = J - 1 i. n! . J - i denotes a multinomial coefficient, j = (J\, . . . ,j n ) T G Nq, and 



{n—k+l n—k+l 
jeN n-k+l. £ iji=nand ji = k\. (7) 

i=l i=l 



Alternatively, one can use Bell polynomials to reformulate (|6]). These are defined by 

(\ n—k+l , v ji 
,1 n . . J n (?) • (8) 
31, • • • ,3n-k+lJ l=1 

This implies that ([6]) can be written as 

n 

(/ o g)W(x) = f {k \ 9 {x))B n ^g'{x),g"{x), g^ k+1 \x)) (9) 
k=l 

In the sections to come, we frequently need the following results. Here, (x) n = 
x ■ {x — 1) • . . . • (x — n + 1) denotes the falling factorial, and s{n, k) and S(n, k) denote 
the Stirling numbers of the first and second kind, respectively, given by the recurrence 
relations 

s{n + 1, k) = s(n, k — 1) — ns{n, k), 
S(n + 1, k) = S{n, k - 1) + kS{n, k), 

for all k G N, n G N , with s(0,0) = 5(0,0) = 1 and s(n,0) = s(0,n) = 5(n,0) = 
5(0, n) = for all n G N. Note that for n G N (in particular n^0), the Stirling numbers 
of the first kind satisfy 

n 

{x) n = Y J s{n,j) X i. (10) 

3=1 

Lemma 3.1 

Let -B n] fc be the Bell polynomial as in {|Sp and n G N. Then 



n; 



(1) For j G 7> n , fc , ELi (* " Oil = 

(2) B n:k (x, ...,x) = S(n, k)x k , k G {0, . . . , n}; 

(3) £ n , fe (-x, . . . , (-l)"- fc+1 x) = (-l) n S(n, k)x k , k G {0, . . . , n}; 

(4) sign(B n ^(g' (x) , g" (x) , . . . , g( n ~ k+1 > (x))) = (— \) n ~ k for all x if g' is completely mono- 
tone. 
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Proposition 3.2 

Let 



(1) For all k € {0, . . . ,n}, B n , k ((x) x y x -\ (x) n . k+1 y x -^~ k +^) = y xk - n s nk (x); 

(2) ELi("lW(x) = {-x) n - 

(3) If x € (0, 1], sign( Snfc (x)) = {-l) n -K 



3.3 The main result 

We are now able to derive a general formula for the derivatives of the inner generators 
and also for the density of nested Archimedean copulas of Type Q . It will follow from 
Faa di Bruno's formula that the derivatives of the inner generators ipo s (t\ Vq) are the inner 
generators themselves times a polynomial in — vq. The product of these derivatives can 
then be computed as a Cauchy product. Interpreting the appearing quantities correctly 
allows us to compute the expectation with respect to Fq via the derivatives of ipo . This 
solves all three of the above challenges. 

Theorem 3.3 

Let i^s €. ^oo, s € {0, . . . , do}, such that ifio s has completely monotone derivative for all 



s(n, l)S{l, k)x l = (-l) n 53|s(n, l)\S(l, k){-x) 1 . 



l=k l=k 



Then 



s e {1,.. .,d }. 
(1) For all n G N, 



n 



(11) 



k=l 



where 



with sign(a Sinfc (t)) = (-l) n k 
t. 

(2) The density of ([I]) is given by 




1 for all 



(12) 



= e b d d ° k (t(u))4 k \t( U ))) . n n^^'K-) 



(13) 




where 



t(u) = (ti(tti),. . .,tdo(«do)) 
do 



#*(*(«)) = E U^dsjsiUUs)) 



(14) 
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with d = (di, . . . , dd ) T and 



Q§k = \j e Ndo : I> = fc > ^ <d„ s€{l,..., d }}; 

that is, 6^ is a coefficient in the Cauchy product of the polynomials J2t=i a s,d s k(t)(—vo) 
Remark 3.4 



1) We see from Theorem |3 . 3| Part (1) that all derivatives of ipo s (t; vq) are of similar form 



in vq, namely ipo s (t; vq) times a polynomial in — vq where the coefficients a s ,d s fc(t s (w s )) 
are the Bell polynomials evaluated at the derivatives of the nodes ipos- This structure 
is crucial for solving Challenge 3 since it allows one to compute the expectation with 
respect to Fq explicitly. 

(2) We see from Theorem 3.3 Part |(2)| how the (log-)density can in general be evaluated. 
It involves th e sign-a djusted derivatives of ipo which are known in many cases; see 
Hofert et al. (2012b). Furthermore, the quantities 6^° fc , k £ {do, . . . , d}, have to be 



computed. The remaining parts are comparably trivial to obtain. 
(3) If there are degenerate child copulas, that is, there exists a subset S of indices such 



that d s = 1 for all s £ 5, then a straightforward application of Theorem 3.3 (1) 
shows that 



-d s \ do d s 

en = ( e b%(t(u))^ k+ds \t(u))) . n nwr 1 )^ 

= d' 7 8=1 j=l 



where ds = J2sesds = \<S\, d' = do — ds and d' is the vector containing all the 
dimensions d s for s ^ S. 



4 Example families and transformations 

4.1 Tilted outer power families, Clayton and Gumbel copulas 

In order to construct and sample new nested Archimedean copulas it turns out to be 
useful to consider certain generator transformations; see Hofert (2010) for more details. 
One such transformation leads to tilted outer power generators 



^(t) = V((c 9s +t) 1/9s -c), 
for a generator ip £ ^oo, c 6 [0, oo), 9 S 6 [1, do), and s G {0, 



(15) 

, do}. Note that generators 
of this form are elements of ^qq. It follows from Equation ^ and Proposition 3.2 Part 

. + t) that the derivatives of t/io are 



1) (with x = l/6>o and y 

n 

4 n \t) = E V> (fc) ((c 9 ° + - c)(c e ° + t) k ^- n s nk {l/9o) 



(16) 



k=l 
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For nesting generators of Type (U5I), the nodes are given by 



^Mt) = (<A+t)°" 



Jo 



This implies that tilted outer power generators of Type (15) fulfill the sufficient nesting 
condition if 9q < 6 S . Furthermore, 

^o1 ) (t) = (a s ) k (c e °+tr°- k , ken. 



By Proposition 3.2 Part (1) (with x = a s and y = c 3 + t), this implies that 



\a s k—n , 



By Theorem 3.3 Part |(1)| this implies that the inner generator 

^ 0s (t;v ) = exp{-v ((c 6 ° + t) a ° - c e °)), 

has derivatives 

n n 



47(t;v ) = iPos(t;v ) E asM^-vo)" = *M*;«o) E^ + t) a ° k - n s nk (a s )(-v ) k . 
k=l k=l 

(17) 



Note that (^J 1 )'^) = O s (ip- 1 )'(u)(c + if)' 1 (it)) 93 ' 1 . By Equation ([16]), Theorem 
Part |(2)| and slight simplifications, we thus obtain 

<<«) = ( E &&(*(«)) (Z^^ + t(u)) 1 ^ -c)(c e ° +t(u)y^- k s kj (l/9 



3.3 



k=do 

do d s 

8=1 j=l 



(18) 



Remark 4.1 (Clayton and Gumbel copulas) 

(1) By taking ip(t) = 1/(1 + t) and c = 1 we see that the tilted outer power generator 
(15) is ip s (t) = (1 + t)^ 1 ^ 3 , that is, a generator of the Clayton family. As a special 
case of this section, we thus obtain the inner generator derivatives and the densities 
of nested Clayton copulas. Concerning the former, we obtain from (17) that 

n 
k=l 



Concerning the latter, plugging in the corresponding quantities in ( 18 ) and simplifying 
the terms (in particular, the power of 1 + iI)q 1 {C{u)) can be taken out of the inner 
sum), we obtain 

= ( e + *wr (fc+1/0o) E(-i) fc -''Mi/0o) 

V k=d j=l 
d , d a N -(l+0„) 



n# n 



3=1 
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By Proposition 3.2 Part |(3)[ we can further simplify this expression and obtain 



c(u)=(-i) d ( x: 

^ k=d 

do / d a x -(i+e a ) 



n# n 



'■SJ 



for the density of nes ted Clay ton copulas of Type Q. This formula also follows 



directly from Theorem 



3.3 



Part 



(2) 



-l/0 o )fc(l + t)- (fc+1/e ° r and simp 



by plugging in the generator derivatives ipQ (t) = 
ifying the expressions. 

(2) Interestingly, also the inner generator derivatives and the densities of nested Gumbel 
copulas of Type Q follow as special case of nested tilted outer power families. To 
see this take ip(t) = exp(— t) (the generator of the independence copula) and consider 



a zero tilt (so c = 0). It follows from (17) that 

n 

4"\t;v ) = ^os(t;v ) J2 s nk (a s )t a ° k - n (-v ) k . 

k=l 

Concerning the density, a short calculation shows that 

c{u) = { - 1)d fM( S ^))(EH tt ) v %(^ 

\ ' k=dn 7=1 



do / d s 

U 0d s s n- io s'' 



8=1 



■V 



s-1 



where C is ([I]) and II is the independence copula (hence the product of its arguments). 
As before, this result can also be directly obtained from 3.3 Part 1(2)] based on Gumbel's 



generator derivatives il^Q k \t) = (ipo(t) / t k ) J2j=i s kj(^/^o){— ^/^oV as derived in Hofert 
et al. ||2012bb. 



4.2 Ali-Mikhail-Haq copulas 

A nested Archimedean copula of Type ([I]) with all components C s , s £ {0, ...,drj}, 
belonging to the Ali-Mikhail-Haq family is a valid copula according to the sufficient 
nesting condition if 8q < 6 S for all s S {1, . . . , do}. The generator ipo s {t] vq) is given by 



1 



l-0o)exp(i) 



Bo) 



I'D 



\exp(t) - 9 0s 



VQ 



where 6q s 



0q)/(1 — 9q) 6 [0, 1) and vq G N. It can be shown from (|9j), Lemma 3.1 



(2)1 and (gOj) that 

n 

4?(*5«o) = i>0s{t;v ) «nfc(l/(l - ^o s exp(-t)))(-7j ) fc , 



(19) 



fc=i 
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which reveals that a s>n k(t) = s n k(l/(l — 9o s exp(—t))) in (11). 
Hofert et al. (|2012b|) showed that 



4 k \t) = (-l) fc ^^Li_ fc (0 o exp(-t)), t G (0,oo), k G No, 



1 



where Li s (z) denotes the poly logarithm of order s at z. It follows from Theorem 3.3 Part 
p)] that 



c(u) = (-1) 



d 1 ~ #o 



E b d d y(t(u))(-l) k U. k (e exp(-t(u))) 



k=do 



do / d s \ 

na-a-n nMi-^(i-«.i))j 

s=l V i=l y 



4.3 Joe copula 

The inner generator and its derivatives 

Nested Joe copulas of Type are valid copulas if 6q < @s for all s G {1, . . . , do}. The 
generator ipo s (t;vo) is given by 



^ te (t;«o) = (l-(l-exp(-t)) a 'r, 



(20) 



where a s = 0q/9 s , s G {1, ... , do}, and -uo G N. 
A rather lengthy calculations shows that 



m 



4:\t;v ) = i> 0s (t;v )(-l) n J2 S(n, 

m=l 

n 

= lp0s{t] Vq) ^ CLs,nk(t)(-V Q ) k 

for 

asMV = ( ("I)""* E 5 ( n > 



t \ m m m 



1 - e~* 



k=l l=k 



+ X 



k=l 



, m 



m=k 



exp(-t) 



1 — exp(— t) 



^2s(l,k)s m i(a s 



l=k 



where ipy as (t) = 1 — (1 — exp(— t))° s denotes Joe's generator with parameter l/a s = 9 s /6q. 
This is precisely the form as given in (11). 
It follows from Hofert et al. p012bl that 



.(1 -exp(-t)) 1 /^ 



exp(— t) 
1 — exp(— t) 



, t G (0,oo), n G N, 
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where P^ e (x) = Ya=i S(k, /)(/ — ! — 1/9q)i-ix 1 . We obtain from Theorem 



that 



3.3 



Part 



(2) 



c(u) = ^(l -ex P HH))W £ (- 1)fe ^( tH )P^( «*('(«» 

d o ds CI _ .\0 S -1 
ad s TT V 1 u sj> 



exp(— t(u)) 



Note that exp(-i(u)) = IT* i(l - ( x - 



4.4 Frank Copula 

Nested Frank copulas of Type ([I]) are valid copulas according to the sufficient nesting 
condition if 9q < 9 S for all s G {1, . . . , do}. The generator i/jo s (t; vo) is given by 

•l_(l_p s e Xp (-£))^° 



Po 



where a s = 9q/9 s , pj = 1 — e 9j , j G {0, s}, and vo G N. Note that this inner generator 
is a shifted (and appropriately scaled) inner Joe generator, that is, 



V'cUMo) 



where h = —logp s ; see Hofert (2010 p. 104) for more details about such generators. 
In particular, with the representation for the generator derivatives for the inner Joe 
generator, this implies that 

,/■(") r/- ? - ^ - ^cL^ + ^o) _ *l>b(h + t;v ) A j fc 

^0s (*) ^Oj - Tj-TT s - Tj-77 2^ a s,nfcl* + "<){-Vq) 



k=l 



fp0s(t; v ) a l,nk(t + h)(-v ) 



k=l 

and thus that a s>n fc(t) = al nk (t + h), that is, the coefficients of the polynomial in —vq 
for the derivatives of the inner Frank generator are the ones of the inner Joe generator, 
appropriately shifted. 

It follows from Hofert et al. (|2012b) that 

4 k \t) = (-l^LL^poexpH)), t G (0,oo), k G N . 



Theorem 3.3 Part |(2)|then implies that 



c(u) = (-l) d ( £ ft(t(«))(-l) fc Li_ (fc _ 1) (p exp(-tH)) 
^k=d 

exp(-9 s u sj ) 



do c 
s=l 3=1 X 



exp(-9 s u 



S3 ) 
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4.5 A nested Ali-Mikhail-Haq o Clayton copula 

If -00 is the generator of an Ali-Mikhail-Haq copula and ip s , s £ {1, ... , do}, generate 
Clayton copulas, then Hofert ( 2010[ p. 115) showed that the sufficient nesting condition 
holds if 6 a S [l,oo), s G {1, . . . ,do}, so one can build nested Archimedean copulas of 
Type ([I]) with the root copula Co being of Ali-Mikhail-Haq and the child copulas C s 
being of Clayton type under this condition (referred to as Ali-Mikhail-Haq o Clayton 
copulas). In this short calculation shows that 

iP 0s (t;v Q ) = i>{{l + t) 1 ^ -1) 



for ip(t) = (1 + (1 — 6o)t) v °. We can thus apply (16) with c = 1 to see that 

n 

4l\t-,v ) = E^ (i) ((! + t ) 1/ds - m+ty /es - n s nj (i/9 s ), 

3=1 

where 

^•)((i + t )V». - 1) = (1 - 9 y^ 0s (t;voHos(t;j) E s(j,k)(-v ) k . 



(21) 



k=l 



Plugging this result into (21) and interchanging the order of the two summations, we 
obtain 

n , n s 

4:\t; v ) = ^ 0s (t; «o) E E k)s nj (l/e s )^ 0s (t; - 6 Q y(l + t^ e '~ n (-v ) k , 
k=i^j=k ' 



which provides the structure of the coefficients a s ^ n k in (11 ), namely, 

3 



a s , nk {t) = E ] s(j, k)s nj (l/9 S/ yeo+ + 

j—k 



1 



70 



(1 + *) 



j/Vs-n 



It is clear from Theorem 3.3 Part |(2)| that the density for the nested Ali-Mikhail-Haq 
o Clayton copula basically consists of the corresponding pieces of the Ali-Mikhail-Haq 
and the Clayton density we have already seen earlier. It is given by 



c(u) = (-l) di ^( E bi k (t(u))(-l) k U_ k (9 eM-t(u)))) IK 



k=do 



do 



8=1 



3=1 



-(i+e s ) 



d ° C (u 

Note that exp(t(u)) = 7 — 



8=1 



1-C s (u s )) 
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5 Numerical evaluation 
5.1 The log-density 

In statistical applications one typically aims at computing the log-density. From a 
numerical point of view, this is typically not as trivial as computing the density and 
taking the logarithm afterwards. Often, the density can not be computed without running 
into numerical problems, hence taking the logarithm of the density faces the same problem. 
However, an intelligent implementation of the log-density is possible (and often even 
required), see the implementation in the R package copula. 

We now briefly explain how one can efficiently compute the log-density of a nested 
Archimedean copula of Type ([I]). Recall from (13) that 

j d s do d a 

= e b% k (t(u))4 k \ t (u))) . n n^/K-)- 

^ k=d ' s=l j=l 

Let us first think about the signs of the terms b^ k (t(w)), k £ {do, . . . , d}. By Theorem 



3.3 



(1) we know that sign(a St d s j s (t s (u s ))) = (-l) da js , thus 



sign[[a s4sjs (t s {u s )) = (-l)s=i . 

s=l 

Recall from ( 14| the structure of b d ^ k {t{u)), which is the sum in j G Q'dk over Hs=l a s,d s js 
For such j, it follows from the definition of Q^fc that Es=iis = ^- Furthermore, note 
that J2t=i ds = d, hence 

sign6j fc (tH) = (-l) d - fc . 

This implies that 

(d do d 3 

e (-i) d - fc ^ fc (tH)(-i)vs fe) (^)) • n nc-v^'K-) (22) 
k=d ' s=l j=l 

where we note that ]T*L i ]lfc s = i(-l) = (-l) d - 

We see from ( 22 ) that all appearing quantities are positive which is quite convenient 
for computing the log-density 

/ d s d d s 

io g cH = iog ^(-i)"^(t( U ))(-i) fc 4 fc) (*(u)) +EE 1 °g(-^- 1 ) , (^). 

^k=d ' s=\j=\ 

Since the latter double sum is typically trivial to compute, let us focus on the first sum. 
To compute the (intelligent) logarithm of this sum, let 

x k = log((-l) d - fc 6* (t(tt))) + log((-l) k 4 h \t(u))) , k G {do, ...,d}, 
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and note that 

log £ (-l) d -%: k (t(u))(-l) k 4 k \t(u)) = log £ exp(x fc ) 
k=do k=do 

d 

= Smax + log eX P( X k ~ %ax), 

k=dg 

where x max = max^^Krf Xk- Since all summands in the latter sum are in (0,1], the 
corresponding logarithm can easily be computed. It remains to discuss how the Xk, 
k G {do, . . . , d}, can be computed. 

For computing the Xk, k E {do, . . . ,d}, efficient implementations for the functions 
log((—l) k t/jQ k \t)) in the R package copula can be used. Computing the quantities 



log((— 1) bj ] k (t(u))) is more challenging. Recall from (14) that 



«&(*(«)) 



do 

n 



(**(««))> 



where a S) d a j a {t s (u s )) is given in (12). For computing the function Sd s j s that often 



appears in a St d a j a {t s { u s)), the function coeffG in copula can be used; to be more 
precise, (-1) ~ (d s -j s ) * copula: : : coeffG (d s ,x) computes s<i s j s ( x )- For the summation 
over the set Qj° fc , the R package partitions provides the function blockparts. With 
blockparts(d-rep(lL,(io) ,k-do)+lL one can then obtain a matrix with do rows where 



do 



each column gives one j G QJ 



k- 



5.2 The -log-likelihood of two-parameter nested Gumbel copulas 

In this section, we compute the -log-likelihood (based on a sample of size n = 100) of 
two nested Gumbel copulas with parameters 9q and 9\ such that Kendall's tau equals 
0.25 and 0.5, respectively. In order to be able to provide graphical insights, we focus on 
two-parameter nested copulas of the form 



C(u) = Cq{ui,Ci{u 2 , ■ ■ -,u d )) 



(23) 



where d G {3, 10}. 

Note that we obtain nested Archimedean copulas of Type (23) from ([T]) by artificially 
thinking of u\ as a child copula V'o('0r7 1 ( tt i)) °f dimension 1, that degenerate 
child copula. For such s, note that a s ^ s j a {t s {u s )) = a s ,ii{t s (u s )) and t s (u s ) equals i/jq 1 



at the corresponding (one-dimensional) argument, which is u\ in (23). It follows from 
the last statement in Theorem 3.3 (1) that a St d a j s (t s ( u s)) = 1 f° r degenerate children. 
This implies that these terms drop out of the product in (14). The set Q^° fc shrinks 
accordingly since 1 < j s < d s = 1 for degenerate children s. 

Figure [T] displays the -log-likelihoods as level plots for the nested Gumbel copulas as 
described above based on a sample of size n = 100 (note the restriction 



< 



The 
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parameters from which the samples were drawn and the minima based on the grid points 
as displayed in the wireframe plots are included. It is interesting to see the behavior 
of the -log-likelihood in the child parameter 9\ when the dimension of the child copula 
C\ is increased (with Cq and its dimension fixed). As can be seen from the plots, the 
-log-likelihood is easier to minimize in #i-direction. This behavior was already observed 



by Hofert et al. (2012b) for Archimedean copulas and can be expressed by the empirical 



observation that the mean squared error behaves like l/(nd) which is decreasing in d 
for fixed n. To see a similar behavior here for the -log-likelihoods of the nested Gumbel 
copulas is not surprising since the marginal copula for u\ = 1 is the Archimedean copula 
C\. Finally, let us remark that the optimization procedure used to generated Figure [T] 
can also be found in the R package copula as a demo. This can directly be used for 
fitting a two-level nested Archimedean copula to a real life data set. Furthermore, similar 
figures as Figure [T] are provided in the demo for a nested Clayton copula. 



-log-likelihood of a nested Gumbel copula 

+ (e„.8,) T x (e^eO 1 



-log-likelihood of a nested Gumbel copula 

+ (e„,e,) T x (e„, n .e,, n ) T 





C(u) = C (u,, C,(u 2 , u 3 )) n = 100 t(6 ) = 0.25 t(B,) = 0.5 C{u) = C {u i ,C,{u 2 Uio)) n = 100 x(e ) = 0.25 x(e,) = 0.5 

Figure 1 Wireframe plots of the -log-likelihood of a three-dimensional (left) and ten- 
dimensionsl (right) nested Gumbel copula C(u) = Co(ui, Ci(ti2, ^3)) with 
parameters 6q = 4/3 (Kendall's tau equals 0.25) and 0\ = 2 (Kendall's tau 
equals 0.5) based on a sample of size n = 100. 



6 Densities for three- (and higher-) level nested Archimedean 
copulas 



In this section, a density formula analogous to (13) is derived for three-level nested 
Archimedean copulas and extensions to higher nesting levels are briefly addressed. 

When working with three or more nesting levels, it turns out to be convenient to 
(slightly) change the notation used in the previous sections. Consider a three-level nested 



15 



6 Densities for three- (and higher-) level nested Archimedean copulas 



Archimedean copula of the form 

C(u) = Ci(Cii(Ciii(uni), . . . , Ciidu(wild dll ))) • • ■ , Cldi(Cldil(l*ldil)3 ■ ■ ■ i Cldid ldl (wididi, 

(24) 

where u SlS2S3 = {u SlS2S3 i, . . . , u slS2S3dsiS2S3 ) T denotes the argument of C SlS2S3 (the copula 
generated by ip SlS2S3 ), d SlS2S3 denotes the dimension of C SlS2 s 3 , and d SlS2 denotes the 
dimension of C SlS2 (the copula generated by tp SlS2 )- Here and in the following, si always 
equals 1, S2 € {1, . . . , d Sl }, and S3 S {1, ... , d SlS2 }. Note that it is convenient to think of 
(24) as a tree; see Figure |2j Furthermore, let 

de l s 2 s 3 
S4=l 

u s 1 s 2 = { u Sl s 2 li ■ ■ ■ ' U s 1 s 2 d aiS2 ) ) 
^sis 2 ( u siS2) = {tsiS2l{ u sis 2 l)i ■ ■ ■ 1 tsiS2d aia2 ( u siS2d aia2 )) j 
Csis 2 (' Us i s 2) = Cs 1 s 2 (C Sl s 2 l(u Sl s 2 i) , . . . , C SlS2 rf siS2 (w SlS2 rf SiS2 )) 

= V^]^ ( V'siS2,SiS2S3 (^SlS2«3(' U SlS2S3)) ) ' 

V S 3 = 1 y 

a!s 1 s2 

ts 1 s 2 ( u sis 2 ) = '4's 1 s 2 ((~' SlS2 {u SlS2 )) = ?/' SlS2jSlS2S 3 (t SlS2S3 (u SlS2S . i )) , 

^3 = 1 



where 



i J s 1 s 2 ,s 1 s 2 s 3 — 4'sis 2 ° '0S1S2S3 



and C* s denotes the (marginal) nested Archimedean copula with root C SlS2 . Note 
that the dimension of the root copula C SlS2 of the nested Archimedean copula C* 1S2 

is d slS2 which is in general not equal to the dimension d SlS2 . = 2^3=1 d slS2 s 3 of Cs lS2 - 
Furthermore, the root copula C\ (= C Sl ) of the nested Archimedean copula C has 
d% (= d Sl ) arguments, the S2th of which has d SlS2 -many arguments. Overall, d si .. = 
Ss 2 Li Ss!=i d slS2S3 equals d, the dimension of C. 



In order to compute the density c of C, we use a similar idea as in Section 3.1 
By replacing ifti( = V>si) = £<S[i<i] with the corresponding integral, we obtain 

/>00 ^1 

C(u) = / J] exp(-v^-\c: iS2 (u SlS2 ))) dFi(«i) 

J0 S 2 =l 
POO ^1 

= / II ^^^(dsjKiss);^)^^) 

J0 8 2 = 1 

where 

ip Sl , SlS2 {t; vi) = exp(-viip Slt s lS2 (t)) 
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uiidni Miidnd lldll u ldlll u ldl id ldll u ldldldi i ui dldldidldiiid ^ 



J V J 



Figure 2 Tree structure (some arguments are omitted) for a three-level nested Archime- 



dean copula of Type ( 24 ) . 



and thus 



c(«) = / H q- ^.-iwC^aaCM-i-a);^)^^). (25) 

JO S2 = 1 OU SlS2 

where — denotes the derivative with respect to all components of u SlS2 (which is a 



vector of length d SlS2 .). 

Similar as in Section [3.1| we observe the following key challenges: 

Challenge 1 Find the derivatives in the integrand; 

Challenge 2 Compute their product; 

Challenge 3 Integrate it with respect to the mixture distribution function F\ = CS~ l [tpx] . 

We will first solve Challenge 1 by considering a multivariate version of Faa di Bruno's 
formula. For suitable functions / : M — > M and g : M. n — > M, it follows from Hardy (2006) 
that 

f( g ( x )) = Y / f {k \g(x)) £ II rf — ^aix), x = ( Xl ,...,x n ) T , (26) 
ax k=i v . M=kBe *nieB<toi 

where the last sum extends over all partitions 7r of {1 , . . . , n} with k elements and the last 
product over all blocks B of ir. Observe that, if then the univariate 

Faa di Bruno's formula (|9[) can be restated as 



k=l it:\n\=k SStt 

where 7r is a partition of {1, . . . , n}. Comparing this identity with ^ yields 

B n ^g\x),...,g^- k+1 \x))= ]T \{g m \x). (27) 

7r:|7r|=fc Bdir 



This will be used in the following lemma, which is a special case of (26) with stronger 
assumptions on the function g. It will then lead us to a solution for Challenge 1 by 
choosing suitable functions / and g. 
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Lemma 6.1 

Suppose there exists a partition {Bi, . . . , B m } of {l,...,n} with \B{\ = d\ for I G 
{l,...,m} (with J2h=idi = n), such that for any indices k\ G Bi and k 2 £ Bj, for 
i,j G {1, . . . , m} with i ^ j, the partial derivative of g(x\, . . . , x n ) with respect to x kl 
and Xk 2 equals zero, that is 

d 2 

- — - — g(x)=0, for all he Bi, k 2 eBj, i,j e{l,...,m}, i^j. (28) 

Moreover, suppose that for any I G {l,...,m} and any subset B of Bi, there exist 
functions h\\ and h\ 2 such that 

TT ^-<K*) = 4 B|) (M*)) II ^M*)- ( 29 ) 

Then one has, for any suitable function /, 



fc=1 ieQ? fc l=i 



where Q^k ^ s defined as in Theorem 3.3 (2) and d = (d\, . . . , d m ) 



We are now in the position to solve Challenge 1. By applying Lemma 6.1 with 
x = u slS2 , n = d slS2 ., m = d slS2 , B t = {sis 2 ll, . . . , sis 2 ld SlS2 i} (slightly abusing the 
notation), d SlS2 = (d SlS2l , . . . , d SlS2dsiS2 ) T , f(t) = t/> Sl , aia2 (t; v{), g(x) = t* SlS2 (u SlS2 ), 

hn{t) = ip sl s 2 ,s 1S2 i(t), and h l2 (u) = t SlS2 i(u SlS2l ) = E S 4=f ^w^i^sJ) the derivative 
in the integrand in (25) is given by 

d 

"4 } s 1 ,s 1 s 2 {'tsis 2 ( u sis 2 



dllsis 2 



n n (^sXsjius^sssi)) e ^£s 1S2 (*s 1S2 (>s 1S2 );ui) 

8 3 =1 s 4 = l ' 1=1 



d ai s 2 



E IT ^ d s 1 a 2 a z ,js 3 (i'4 , S 1 l 2 ,S 1 S 2 S 3 (ts 1 S 2 S S {u SlS2 s s 

)))*:{ 
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By applying Theorem 3.3 (1)| this derivative can be written as 
d 



du SlS2 

n II (^s 1 1 s 2 s 3 y(Us 1 s 2 s s s 4 )j^s 1 ,s 2 s 2 (t* SlS2 {u SlS2 );V 1 ) 
S3 = l S4=l 
dsis 2 - I 

E J2 a sis 2 M t *sis 2 { u sis 2 ))(-Vl) k 
1=1 k=l 

d si s 2 . 

e n B dsiS2S3)jS3 mii 2 

,s 1 S2S3i^siS2S-i( u sis 2 s: j )))ke{l,...,ds 1 s 2 s a -js 3 +l}) j 

dsiS2 .3=1 

Interchanging the order of summation in the double sum yields 
d 

(JU SlS2 

ds 1 s 2 d si s 2 s 3 

n n (^s^sj ( u sis 2 s 3S4 ) w (w aia2 );«i) 

S3=l S4=l 

E (-^l) fc E (^W^l^Kl^)) 



k=l l=k 

B l s 2 



U/B l s 2 \ \ 

,SlS 2 S 3 (^lS 2 S3(' U 'SlS 2 S3)))fce{l,...,d i , 1 s 2 S3-is3+l}) J J 

-,<*si s 2 S3=l 



which can be written as 



d 

«- — ^«i,-i« 2 (<; i «(w«i- 2 );«i) 



is 1 s 2 rfs 1 s 2 s 3 



(u slS2 );vi) 

i3 = l S4 = l / 

! Sl s 2 . 

E a ««,d.i. a .k( t W u <*«a))(-' l, i)*> (30) 
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for 



si «2 s 2 ■ (^si s 2 ('"si s 2 ) ) 



l=k 



E ( a siS2,lk( t *s 1 s 2 ( U siS2)) 
d Sl s 2 

E IT ^ d s 1 s 2 s 3 ,js 3 ((' t Ps 1 S2,s 1 S2S 3 (^s 1 S2S 3 (Us 1 S2S S )))k£{l,...,d SlS2S:i -j 133 +l}) 



With this notation, the connection with Theorem 3.3||(1) is clearly visible. 
Finally, in order to solve Challenge 2 and 3, one introduces 



(31) 



b d s \ ,k (*(")) = E II a siS2,ds 1S 2-js 2 ( t SlS2( U S 1 S 2 )), 



(32) 



with 



t(u) 


= (*i(wi) T ,-. 




= (t Sll (u Sll ) J 


d si 





\T\T 



! • • • > ^sid si ( M sia! S:L ) ) j 
T \T 



Similarly to Theorem 3.3 (2) one then obtains (with t(u) = ip 1 1 (C(u)) as before) 



oo i 



d 



c ( u ) = I n jt 1 — ^^^(^Ksjsvij^i^i) 

JO S2= l ou sis 2 



S1 u. Sl s 2 "s lS2 S3 



n n n (^u/k 

82=1 »3=1 S4=l 



E &£,*(*(«)) 



k=d a 



I ( n ^ 1 , slS2 (t: 1S2 K 1S2 );^i))(-^i) fc di ? i(^i) 



n n n (^u)'^-.--*)) e ^fc,*^))^^))- ( 33 ) 



82 = 1 S3 = l 84=1 



Remark 6.2 

The pattern to compute the density of nested Archimedean copulas with more than three 
levels can be deduced from the previous computations, the following heuristic argument 
shows how. In order to understand the reasoning, it is useful to remind ourselves that 
the structure of nested Archimedean copulas can be depicted by trees; see Figure [2] 
for a tree representation of ( 24 ) . Let L denote the number of levels (with ( 24 ) having 



20 



A Appendix 



L = 3). Thanks to our notation, we can easily identify a certain branch of the tree 
with the corresponding sequence of indices. Each time a nesting level is added, for 
each branch s\S2 ■ ■ ■ si, I £ {1, . . . , L}, finite sequences of coefficients a's and 6's (similar 



to a 



siS2,d s 



and b£ k above) will appear and their structure can be deduced from 



Equation (31). More precisely, as in Equations (14) and (32), the sequence of 6's can 



always be interpreted as the coefficients of the Cauchy product of the polynomials with 
a's as coefficients. 

The structure of the a's at each branch s\ . . . si is more complicated. For any branch 
s\ . . . sl, that is on the ultimate level of nesting, the a's are simply the Bell polynomials 
applied to the function tp Sl ... SL _ ltSl --s L an d its derivatives. For any other branch s± . . . si, 
I & {1, . . . , L — 1}, the coefficients a's are the (Euclidean) inner product of the vector 
of all Bell polynomials applied to the function tp S i...si,si...si +1 and its derivatives with 
the vector of all coefficients 6's, the exact structure of Equation (31 ). In Equation (31 ), 
the level I is equal to 1 = L — 2 and the term a SlS2 jk(t* SlS2 (u slS2 )) stands for the Bell 
polynomial applied to VW-»si,si..-*i+i and its derivatives, while the Zth member of the 
sequence of 6's is defined by 



j^ri, 33=1 



II ^d slS2 B 3 ,js 3 {(^s 1 l 2 ,s 1 S2S 3 ( t s 1 S2S 3 {u SlS2 s 3 ))) k ^ 1 ^ ^ d ^^ 



-js,+l} 



the term appearing in (31) 



Acknowledgements 

The authors would like to thank Paul Embrechts (ETH Zurich) for giving us feedback 
and the opportunity to work on this topic. 



A Appendix 

Proof of Lemma 13.11 

(1) Er=i fe+ v - u = x Er=i fc+1 31 - Er=t +1 & = xk- n . 

(2) The identity B n ^{l^ . . . , 1) = S(n, k) can be found, for example, in Comtet (1974 
135) or Charalambides ( |2005 p. 87). It then follows that 



P- 



B„, t (*,...,*)= V ' 1| 

= x k B ntk (l,...,l) = S(n,k)x k 



n-k+l / x \Jl 



since Y^=\ +l Jl = khy definition of V, 



n,k- 
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(3) By definition of V n ,k it follows from Er=i fc+1 Iji = n and|(2) 
B n , k (-x,...,(-l) n - k+1 x)= E 



that 



\Jli ■ ■ ■ i Jn— k+1 . 



= (-l) n B n , k (x, ...,x) = (-l) n S(n, k)x k . 

(4) By (§, 

^n, fc (5 / (x),5 ,/ (x),...,^- fe+1 )(x))= J2 



n 



\ j\i • • • i Jn—k+l , 



i=i 



l\ 



(34) 



Note that sign(g(')(x)) = (—1)' 1 for all x and I G {1, . . . , n — k + 1}, so that the 
sign of the Zth factor in the product in (34) is (— This implies that 



n-k+l 

sign Yl 
i=i 



g( l \x) 
11 



n— fc+1 n — fc+ 1 n— fc+1 

i« E E lii- E A 

= (-1) 1=1 = (-1) 1=1 !=1 . 



Now since we sum over jf € T^fc, we see from ([7]) that 

' g">(x 



sign ]J 



n-k+l , J$)l x \ \ .11 



1) 



n—k 



and thus the whole sum in (34) has this sign. 



□ 



Proof of Proposition |3.2| 

(1) Let h(t) = exp((yt) x ) . It follows from Faa di Bruno's formula (([6]) and ^ with 
f(x) = exp(x) and g(t) = (yt) x ) and Lemma 3.1 Part |(1)| that 



h^(t)=h( t )Y: e 



n 



\Jl j • • • ; Jn— k+1 



n 

i=i 



(x)w x t 0[ 



n 



,Jl j • • ■ j Jn—k+l 



II 



\Jl > ■ • • ! Jn—k+l 



n-k+l ,( n ,\.fic\ h 



n 



\ n— fc+1 

n 



k=ijeV n! k 

n 

*(*)£ E 

fc=l jGP„, ft 
n 

^)Ey" E 

fc=l jeV nik 

^ E y^fcCCx)!^- 1 , . . . , (x) n _ fc+l2/ -("- fc+1 ))t 
r fc=i 



xk—n 



1=1 



(x)iV 



t 



x—l \ Ji 



f 



xk—n 
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On the other hand, we may differentiate the series ex pansio n of h and use the identity 
e~ x Efclo k l x k /k\ = ELo S(h k)x k , see Boyadzhiev ( |2009| ), with x being (ytf . With 
Tol it then follows that 



oo / n 



00 j-xk—n 1 

h (n) (t) = Y,(xk) n —y* k = - E ( E *( n > W 



fc=0 



fc=0 x i=l 
\xk 



(yt) 



xk 



k\ 



i E 0^ f>'^rp = 6XI) t t)!) E -(», 0^ E 5(1, k)(ytr" 
1 1=1 k=o K - 1=1 k=0 

hit) 

t 



ir £ (V £ 05(z, k)\ t xk = ^ E £ ^ s (™> O^a, *o) t xk 

1 — n V /— u ' " l — 1 V J— u ' 



k=0 l=k 7 k=l v l=k 

Comparing the two representations for leads to the result as stated. 
(2) Since Yjk=i( x )kS(j, k) = x J and by (10), one obtains that 

n n j n 

£(-l)jfcS n fc(a?) = ^s(n,j)x : > ^(-l) k S(j,k) = ^s(n,j)(-xy = (-x) n . 

k=l j=l k=l j=l 



(3) For all x S (0,1], s n k(x) = (—l) n ~ k p(n;k)n\/k\, where the probability mass func- 
tion p(n; k) > (in n £ {k, k + 1, . . . }) corresponds to the distribution function 
whose Laplace-Stieltjes transform is the inner generator appearing in a nested Joe 
copula; consider Hofert (2010 p. 99) with Vq = k and 9o/di = x to see this. This 
representation implies that sign(s n fc(x)) = (— l) n ~ k . 

□ 

Proof of Theorem 13.31 

(1) Apply Q with f(x) = exp(—vox) and g(t) = tpQ S (t). For the statement about the 
signs, apply Lemma 3.1 (4)| For the last statement, note that by Lemma 3.1 (2) 
o a ,n(t) = Bi,i($,.(t)) = = 1) • l 1 = 1 if V s = ^0. 

(2) Given the form ( |11[ ) we see that the product appearing as integrand in ^ can be 
computed via 

do do d 3 do 

II ^0s s) (^(«s); v ) = II E a s,dA t s{u s ))(-v ) k ■ Jl ip 0s (t s (u s ); v ) 



8=1 



=1 k=l 



do 



E bi k (t(u))(-v ) k ■Y[^s(t s (u s );v ). 

k=do 



s=l 



Now note that risLi V'Os^s^s); vq) = exp(— vot(u)). Hence, we obtain 

do d 

H4 d /\t s (u s );v )= E b% k (t(u))(-v ) k exp(-v t(u)). 

s=l k=do 
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By replacing vq by V$ and taking the expectation, one obtains 

do -, d 



E 



H4 d ;\t s (u s );V ) = E b% k (t(u))E[(-V ) k eM-V t(u))] 

k=do 

= E bi k (t(u))4 k \t(u)), 

k=do 



so that, by the density is of the form as stated. 
Proof of Lemma 16.11 

Due to Equation (|28|), observe that the second sum in (26) may be rewritten as 



□ 



e nnSs*-)- e fi( e n 



7r:|7r|=fc -BS7T 



Q\B\ 



Ui£B 9xi 



where tti is a partition of B[ with elements, I G {1, . . . , to}; note that < d;, because 
\Bi\ = d\. Both sides are equal because the remaining terms in the sum on the left-hand 
side are zero as derivatives with respect to two or more variables belonging to different 
partitions vanish. Combining this result with (29), it follows from (26) that 

d 



dx 



k=l 



m*)) = e e n ( e ( n ( h^Hhnix)) n 







da 



Observe that for any partition 717 of Bi, I G {1, . . . , m}, we have that 

n n £m*) = n £rM*)> 



ieBi 



so that 



d_ 
dx 



/w)=E/%w)En n^w e n4 B|) (M-))) 
nn^w e/ w m e n e n^w*))- 



3eQ^ fc i=l7rj:|7r,|=j, BStt, 



Finally, applying Identity (27) completes the proof. 



□ 
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